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Electron  valley,  a  degree  of  freedom  that  is  analogous  to  spin,  can 
lead  to  novel  topological  phases  in  bilayer  graphene.  A  tunable 
bandgap  can  be  induced  in  bilayer  graphene  by  an  external  electric 
field1-5,  and  such  gapped  bilayer  graphene  is  predicted  to  be  a  topo¬ 
logical  insulating  phase  protected  by  no-valley  mixing  symmetry, 
featuring  quantum  valley  Hall  effects  and  chiral  edge  states6  9. 
Observation  of  such  chiral  edge  states,  however,  is  challenging 
because  inter-valley  scattering  is  induced  by  atomic-scale  defects  at 
real  bilayer  graphene  edges10.  Recent  theoretical  work11-13  has  shown 
that  domain  walls  between  AB-  and  BA-stacked  bilayer  graphene 
can  support  protected  chiral  edge  states  of  quantum  valley  Hall 
insulators.  Here  we  report  an  experimental  observation  of  ballistic 
(that  is,  with  no  scattering  of  electrons)  conducting  channels  at 
bilayer  graphene  domain  walls.  We  employ  near-field  infrared  nano- 
metre-scale  microscopy  (nanoscopy) 14-16  to  image  in  situ  bilayer 
graphene  layer-stacking  domain  walls  on  device  substrates,  and  we 
fabricate  dual-gated  field  effect  transistors  based  on  the  domain 
walls.  Unlike  single-domain  bilayer  graphene,  which  shows  gapped 
insulating  behaviour  under  a  vertical  electrical  field,  bilayer  gra¬ 
phene  domain  walls  feature  one-dimensional  valley-polarized  con¬ 
ducting  channels  with  a  ballistic  length  of  about  400  nanometres  at 
4  kelvin.  Such  topologically  protected  one-dimensional  chiral  states 
at  bilayer  graphene  domain  walls  open  up  opportunities  for  explor¬ 
ing  unique  topological  phases  and  valley  physics  in  graphene. 

Hexagonal  two-dimensional  crystals,  such  as  graphene  and  trans¬ 
ition-metal  dichalcogenides,  exhibit  a  pair  of  degenerate  bands  at  the  K 
and  K'  valleys  in  momentum  space.  The  valley  electrons  are  charac¬ 
terized  by  non-trivial  Berry  curvatures,  which  give  rise  to  anomalous 
quantum  Hall  states  in  graphene17  as  well  as  valley  Hall  effects  (the 
valley  counterpart  of  spin  Hall  effects18,19).  Bilayer  graphene  provides 
an  attractive  platform  from  which  to  explore  topological  valley  phys¬ 
ics20:  a  tunable  semiconductor  bandgap  can  be  induced  in  bilayer 
graphene  using  a  vertical  electrical  field,  and  different  topological 
phases  can  exist  in  gapped  bilayer  graphene6,7,11,21-23.  In  particular, 
quantum  valley  Hall  insulator  states  can  exist  in  bilayer  graphene  (with 
any  edge  termination  other  than  the  armchair  type)  with  a  field- 
induced  bandgap,  where  counterpropagating  chiral  electrons  with 
opposite  valley  index  exist  at  their  boundaries6-9.  Such  one-dimen¬ 
sional  chiral  states  are  topologically  protected  as  long  as  there  are  no 
inter- valley  scatterings  generated  by  atomic- scale  defects. 

Domain  walls  between  AB-  and  BA-stacked  bilayer  graphene  pro¬ 
vide  an  attractive  place  to  realize  one- dimensional  chiral  boundary 
states  of  quantum  valley  Hall  insulators  because  a  smooth  domain  wall 
preserves  the  electron  valley  index  in  bilayer  graphene,  unlike  valley¬ 
mixing  atomic  defects  at  graphene  edges.  The  integrated  Berry  curv¬ 
ature  of  bilayer  graphene  valence  bands  for  each  valley  is  characterized 
by  a  quantized  “valley”  Chern  number11-13,  the  sign  of  which  depends 
both  on  the  direction  of  the  vertical  electrical  field  and  the  layer 


stacking  order.  Martin  et  aln  first  proposed  topological  domain 
walls  between  band-inverted  bilayer  graphene  insulators  generated 
by  opposite  vertical  electrical  fields.  Experimental  realization  of  such 
electric-field  walls,  however,  is  extremely  challenging.  An  alternative 
approach  is  to  exploit  different  stacking  orders.  AB-  and  BA-stacked 
bilayer  graphene  under  the  same  vertical  electrical  field  represent  two 
distinct  topological  phases  with  opposite  valley  Chern  numbers.  The 
stacking  boundary  between  AB  and  BA  domains  forms  a  layer  stacking 
domain  wall  (Fig.  la).  Valley-polarized  chiral  electrons  are  pre¬ 
dicted12,13  to  flow  along  such  a  domain  wall  in  gapped  bilayer  graphene, 
where  K  and  K'  valley  electrons  propagate  in  opposite  directions. 
Recent  transmission  electron  microscopy  studies  show  that  smooth 
AB-BA  domain  walls  with  a  width  of  several  nanometres  exist  naturally 
in  suspended  bilayer  graphene24,25,  but  such  suspended  transmission 
electron  microscopy  bilayer  graphene  samples  on  through-holes  (that 
is,  holes  that  penetrate  through  a  substrate)  cannot  easily  be  com¬ 
bined  with  bandgap  control  and  transport  measurements.  Here  we 
use  a  new  approach  to  image  in  situ  AB-BA  domain  walls  in  exfo¬ 
liated  (that  is,  prepared  from  graphite  using  tape)  bilayer  graphene 
on  Si02/Si  substrates,  and  demonstrate  topologically  protected 
charge  and  valley  transport  along  such  boundaries  for  the  first  time. 

Few-layer  graphene  was  mechanically  exfoliated  on  280-nm-thick 
Si02/Si  substrates,  and  the  layer  thickness  is  determined  using  optical 
contrast  measurements.  To  identify  AB-BA  layer  stacking  domain  walls 
in  bilayer  graphene  we  used  near-field  infrared  nanoscopy14,  as  illu¬ 
strated  in  Fig.  lb.  Infrared  light  at  6. 1  pm  was  focused  onto  the  apex  of  a 
metal- coated  atomic  force  microscope  (AFM)  tip  with  curvature  radius 
r  ~  25  nm,  and  the  scattered  infrared  radiation  was  collected  and  mea¬ 
sured  by  an  HgCdTe  detector  in  the  far  field.  The  tip- enhanced  infrared 
scattering  provides  a  local  probe  of  the  material  infrared  responses  with 
—40  nm  spatial  resolution.  The  AB-  and  BA-stacked  bilayer  graphene, 
being  inversion  symmetric  to  each  other,  have  exactly  the  same  elec¬ 
tronic  bands  and  infrared  absorption.  The  local  band  structure  of  the 
domain  wall,  however,  is  strongly  modified  relative  to  the  bulk  and 
exhibits  different  infrared  contrast.  Figure  lc  and  d  displays  AFM 
topography  and  near-field  infrared  nanoscopy  images  of  the  same 
bilayer  graphene  on  a  Si02/Si  substrate.  The  bilayer  graphene  topo¬ 
graphy  shows  a  step  height  of  —  1  nm,  consistent  with  the  bilayer  thick¬ 
ness.  Within  the  bilayer  graphene  flake,  the  topography  is  featureless 
with  a  height  fluctuation  of  less  than  0.2  nm.  In  the  near- field  infrared 
nanoscopy  image,  however,  distinct  features  appear  in  the  bilayer  gra¬ 
phene:  multiple  bright  lines  extend  across  the  bilayer  graphene  flake 
and  terminate  at  the  edges.  These  bright  lines  in  the  infrared  image 
indicate  a  change  of  local  infrared  conductivity  in  the  bilayer  graphene. 
No  corresponding  features  can  be  discerned  in  the  topography  image 
(Fig.  lc),  which  excludes  the  possibility  of  surface  contaminations  or 
graphene  folding  along  these  lines.  We  attribute  these  local  infrared 
response  changes  to  AB-BA  domain  walls  in  bilayer  graphene. 
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Figure  1  |  Imaging  AB-BA  domain  walls  in  exfoliated  bilayer  graphene. 

a,  Illustration  of  a  domain  wall  (shaded  area)  between  AB-  and  BA-stacked 
bilayer  graphene  domains.  The  yellow  (green)  arrow  indicates  the  one¬ 
dimensional  conducting  channel  at  the  domain  wall  from  K-valley  (K' -valley) 
electrons,  b,  Illustration  of  the  near- field  infrared  nanoscopy  measurement  of 
graphene  on  a  Si02/Si  substrate  (pink/grey).  Infrared  light  (yellow)  at  6.1  pm 
was  focused  onto  the  apex  of  a  metal-coated  AFM  tip  and  the  backscattered 
infrared  radiation  was  collected  and  measured  by  an  HgCdTe  detector  in  the  far 


b 


field,  c,  AFM  topography  map  of  a  bilayer  graphene  sample  on  a  Si02/Si 
substrate,  showing  no  surface  wrinkles  or  defects.  (The  small  triangular  area  at 
the  bottom  right  corner  corresponds  to  a  graphene  monolayer.)  d,  Near-field 
infrared  nanoscopy  image  of  the  same  bilayer  graphene  sample  as  in  c.  Bright 
line  features  across  the  bilayer  graphene  flakes  are  observed  that  are  absent  in 
the  topography  image.  These  bright  lines  arise  from  AB-BA  domain  walls  in 
bilayer  graphene  because  the  domain  walls  have  different  local  electronic 
structure  and  infrared  responses. 


Qualitatively,  we  can  approximate  bilayer  graphene  in  the  domain 
wall  region  as  randomly  stacked  bilayers,  which  has  stronger  local 
infrared  response  than  bulk  AB-stacked  bilayer  graphene  and  leads 
to  the  bright  lines  observed  in  infrared  nanoscopy  images  (see 
Methods  for  details). 

To  confirm  that  the  bright  lines  within  the  bilayer  graphene  are  layer 
stacking  domain  walls,  we  systematically  examined  near-field  infrared 
images  of  exfoliated  graphene  flakes  of  different  layer  thicknesses. 
Figure  2a-c  shows  infrared  nanoscopy  images  of  representative  mono- 
layer,  bilayer,  and  trilayer  graphene  flakes,  respectively.  The  infrared 
images  of  monolayers,  like  that  in  Fig.  2a,  always  show  homogenous 
responses  within  the  flake,  although  bright  features  can  be  observed 
parallel  to  the  edges  owing  to  the  excitation  of  graphene  plasmons.  The 
absence  of  any  structural  domain  boundary  in  monolayers  is  consist¬ 
ent  with  the  high  quality  of  the  exfoliated  samples.  In  bilayer  graphene, 
about  30%  of  the  exfoliated  flakes  show  bright  lines  across  the  interior 
region  of  graphene  like  those  in  Fig.  2b,  and  areas  between  these  lines 
have  the  same  contrast.  We  identify  these  bright  lines  as  AB-BA 
domain  walls  because  they  are  the  only  one- dimensional  structure  that 
exists  in  bilayer  but  not  monolayer  graphene.  This  identification  is 
further  bolstered  by  the  infrared  images  of  graphene  trilayers,  such 
as  that  shown  in  Fig.  2c. 

Bright  lines  extending  through  the  graphene  flakes  are  also 
observed,  but  unlike  graphene  bilayers,  the  domains  separated  by 
the  bright  lines  can  have  different  infrared  contrast.  This  is  because 
more  domain  stacking  orders  can  exist  in  the  trilayer,  such  as  ABA, 
BAB,  or  ABC  stacking.  Although  ABA  and  BAB  domains  are  inversion 
symmetric  to  each  other  and  have  the  same  electronic  structure,  ABC 
domains  have  distinct  band  structures26,  infrared  responses,  and 
Raman  spectra27.  Therefore,  a  domain  wall  between  ABA-  and  ABC- 
stacked  trilayers  will  separate  two  domains  having  different  infrared 
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contrast  in  the  near-field  infrared  image,  as  seen  in  domains  I  and  II  in 
Fig.  2c.  In  addition,  we  performed  micro-Raman  measurements  on 
domains  I  and  II,  and  show  the  corresponding  spectra  in  Fig.  2d.  It  is 
apparent  that  the  two-dimensional  Raman  features  in  domains  I  and  II 
are  different,  and  they  have  been  shown  to  correspond  to  ABC-  and 
ABA-stacked  trilayer  graphene27,  respectively.  It  demonstrates  unam¬ 
biguously  that  the  bright  line  features  in  trilayers  are  layer  stacking 
domain  walls,  and  supports  our  assignment  for  bilayers. 

After  identifying  the  layer  stacking  domain  walls  in  bilayer  graphene 
on  a  Si02/Si  substrate,  we  fabricated  dual-gated  bilayer  field-effect  tran¬ 
sistor  devices4,5  across  the  domain  walls.  As  references,  we  also  fabri¬ 
cated  dual -gated  devices  without  domain  walls  using  the  same  device 
configuration  on  the  same  bilayer  graphene  flakes.  Figure  3a  displays  the 
optical  microscopy  image  of  a  representative  device.  The  white  dashed 
line  outlines  the  bilayer  graphene  flake,  and  the  green  line  denotes  the 
position  of  the  AB-BA  domain  wall  in  the  bilayer  graphene.  Source  and 
drain  electrodes  (electrodes  1-5)  composed  of  40  nm  Au  and  0.5  nm  Cr 
layers  were  defined  by  electron-beam  lithography.  Electrodes  1  and  2  are 
connected  by  a  layer  stacking  domain  wall,  as  are  electrodes  2  and  3. 
Electrodes  3  and  4  and  electrodes  4  and  5,  on  the  other  hand,  are  within  a 
single  bilayer  graphene  domain  and  act  as  reference  devices.  On  top  of 
the  bilayer  graphene  flake,  a  35-nm-thick  A1203  dielectric  layer  was 
grown  using  atomic  layer  deposition.  Last,  the  top  gate  electrodes  of 
width  1  pm  (TGI  and  TG2)  were  patterned  using  electron-beam  litho¬ 
graphy  and  metal  deposition.  Figure  3b  illustrates  the  side  view  of  the 
device  between  electrodes  2  and  4,  with  a  domain  wall  (green  line) 
between  electrodes  2  and  3  and  single- domain  bilayer  graphene  between 
electrodes  3  and  4.  Here  the  doped  silicon  substrate  acts  as  a  bottom  gate. 
By  varying  the  top  gate  voltage  Vtg  and  the  bottom  gate  voltage  Vbg,  we 
can  independently  control  the  electric  field  above  and  below  the  bilayer 
graphene  to  generate  a  tunable  semiconductor  bandgap5. 
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Figure  2  |  Infrared  nanoscopy  images  of  domains  and  domain  walls  in 
exfoliated  mono-,  bi-  and  trilayer  graphene,  a,  Infrared  image  of  a  graphene 
monolayer.  The  darker  region  corresponds  to  the  Si02/Si  substrate. 
Monolayer  graphene  shows  homogeneous  contrast  across  the  whole  sample 
with  no  bright  lines  in  the  interior  area.  (The  bright  lines  parallel  to  the  edge 
arise  from  plasmon  excitations.)  The  same  behaviour  is  observed  in  all 
exfoliated  monolayers,  b,  Infrared  image  of  a  graphene  bilayer,  which  shows 
bright  line  features  across  the  flake.  Such  bright  lines  arise  from  AB-BA 
domain  walls  in  bilayer  graphene  and  they  are  observed  in  —30%  of 
exfoliated  bilayer  graphene  samples.  The  bulk  AB  and  BA  domains  have  the 

Figure  3c  displays  the  electrical  transport  at  4  K  in  a  reference  bilayer 
graphene  device  without  a  domain  wall.  Each  trace  shows  the  gate- 
dependent  resistance  as  a  function  of  Vtg  at  a  fixed  bottom  gate  voltage, 
Vbg.  Vbg  was  varied  in  steps  of  10  V  to  produce  the  different  traces.  The 
resistance  peak  in  each  trace  corresponds  to  a  charge  neutral  point 
(CNP),  at  which  the  Fermi  level  resides  within  the  semiconductor 
bandgap.  CNPs  can  be  realized  with  different  combinations  of  Vtg 
and  Vbg,  but  each  with  a  different  vertical  electrical  field  and  induced 
semiconductor  bandgap  in  bilayer  graphene.  The  zero -bandgap  state  is 
realized  at  Vbg°  =  0  V  and  Vtg°  =  1.8  V,  which  has  the  lowest  peak 
resistance.  The  average  displacement  field5  across  bilayer  graphene  is 
described  by  D  =  V2(Dh+Dt).  Here  Dh  =  £Sio2 1  Vbg  -  Vbg°  |/dSio2,  and 
A  =  £ai2o3  |  Vtg  —  Vtg0 1  / dA\2o3 ,  where  £  and  d  are  the  dielectric  constant 
and  thickness  of  the  respective  oxide  dielectric  layers.  The  bandgap 
increases  monotonically  with  the  average  vertical  displacement  field  D, 
which  results  in  a  larger  peak  resistance  at  the  CNP.  The  experimental 
values  of  the  peak  resistances  are  lower  than  that  expected  from  the 
semiconductor  bandgap  in  ideal  bilayer  graphene,  presumably  owing 
to  impurities  and  defects  in  our  devices.  Nevertheless,  the  resistance 
can  reach  —90  kO  at  4.2  K,  similar  to  that  observed  in  previous  dual¬ 
gated  bilayer  graphene  devices  on  Si02/Si  substrates4. 

Electrical  transport  through  bilayer  graphene  with  an  AB-BA 
domain  wall  exhibits  a  different  behaviour,  as  shown  in  Fig.  3d. 
When  the  bilayer  graphene  bandgap  is  small  (purple  traces),  the 
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same  electronic  structure  and  infrared  contrast,  c,  Infrared  image  of  a 
graphene  trilayer,  which  exhibits  domains  with  different  infrared  contrast  as 
well  as  bright  lines  between  domains.  Different  stacking  orders,  such  as  ABA 
and  ABC  stacking,  can  lead  to  different  electronic  structures  and  infrared 
contrast  in  trilayers,  d,  Raman  spectra  taken  at  domains  I  and  II  in  c;  co  is  the 
Raman  shift.  The  two-dimensional  peak  of  domain  I  (black  curve)  and 
domain  II  (red  curve)  are  characteristic  of  ABC  and  ABA  stacking, 
respectively.  It  confirms  that  the  bright  and  dark  trilayer  regions  in 
c  correspond  to  different  layer  stacking  orders,  and  that  the  bright  line 
features  arise  from  layer  stacking  domain  walls. 

gate-dependent  resistance  curves  are  similar  to  those  in  single-domain 
bilayer  graphene  devices  because  bulk  current  flowing  through  the 
domain  dominates.  With  higher  displacement  fields  and  larger  bilayer 
graphene  bandgaps,  the  peak  resistance  in  each  curve  first  increases 
and  then  saturates  at  — 14  kf}  (red  traces  in  Fig.  3d),  in  striking  contrast 
to  the  monotonically  increasing  resistance  to  over  80  kO  in  the  ref¬ 
erence  device  without  a  domain  wall  (Fig.  3c).  Figure  4a  plots  the 
comparison  of  bilayer  graphene  peak  resistances  in  the  device  with  a 
domain  wall  (red  circles)  and  without  a  domain  wall  (black  squares)  as 
a  function  of  the  vertical  displacement  field  D,  explicitly  showing  that 
saturation  of  the  peak  resistances  occurs  at  a  much  lower  value,  owing 
to  the  presence  of  the  domain  wall. 

These  distinct  transport  behaviours  in  bilayer  graphene  with  and 
without  AB-BA  domain  walls  were  observed  in  all  of  the  bilayer  gra¬ 
phene  devices  that  we  studied  (see  Methods  for  details.)  It  is  apparent 
that  the  conductance  is  always  much  higher  in  bilayer  graphene 
devices  with  AB-BA  domain  walls.  In  addition,  we  performed  temper¬ 
ature  dependence  studies  of  CNP  resistance  in  bilayer  graphene  with 
and  without  domain  walls  from  1.8  K  to  50  K.  We  observed  that  the 
resistance  of  bulk  bilayer  graphene  increases  many  times  at  lower 
temperature,  characteristic  of  an  insulating  phase.  The  resistance 
across  the  domain  wall,  however,  behaves  like  a  metallic  system  with 
much  weaker  temperature  dependence  (see  Methods  for  details).  Our 
results  demonstrate  unambiguously  that  one- dimensional  conducting 
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Figure  3  |  Dual-gated  field-effect  transistor  devices  on  bilayer  graphene 
with  AB-BA  domain  walls  and  their  electrical  transport  properties. 

a,  Optical  micrograph  of  a  dual-gated  bilayer  device.  The  white  dotted  line 
outlines  the  bilayer  graphene  flake,  and  the  green  line  denotes  the  AB-BA 
domain  wall.  Source  and  drain  electrodes  are  labelled  from  1  to  5.  Electrodes  1 
and  2  are  connected  by  a  layer  stacking  domain  wall,  as  are  electrodes  2  and  3. 
Electrodes  3  and  4  and  electrodes  4  and  5,  on  the  other  hand,  are  within  a  single 
bilayer  graphene  domain  and  act  as  reference  devices.  TGI  and  TG2  denote  top 
gate  electrodes,  b,  Side  view  of  the  dual-gated  bilayer  device  between  electrodes 
2  and  4  in  a.  c,  Gate-dependent  resistance  for  bilayer  graphene  without  a 


domain  boundary  (electrodes  3  and  4).  The  backgate  voltage  Vbg  is  varied  from 
0  V  to  140  V  with  a  step  of  10  V  (colours  here  and  in  d  indicate  different  curves), 
and  the  top  gate  voltage  Vtg  is  swept  continuously.  The  resistance  peak  in  each 
trace  corresponds  to  a  CNP.  The  resistance  at  the  CNP  increases  continuously 
(dashed  grey  line)  to  over  80  kQ  at  high  Wg  owing  to  the  opening  of  a  bandgap. 
d,  Gate-dependent  resistance  across  bilayer  graphene  with  a  domain  boundary 
(electrodes  2  and  3).  The  resistance  at  the  CNP  first  increases  and  then  saturates 
at  ~14kQ,  although  the  bilayer  graphene  bandgap  keeps  increasing  with 
increasing  Vhg.  This  much  lower  resistance  signifies  the  presence  of  a  one¬ 
dimensional  conducting  channel  at  the  AB-BA  domain  wall. 


channels  are  present  at  the  AB-BA  domain  walls,  and  that  they  dom¬ 
inate  the  electrical  transport  when  bulk  bilayer  graphene  domains 
become  insulating. 

To  examine  the  conducting  channels  through  the  bilayer  graphene 
domain  walls,  we  systematically  investigated  bilayer  graphene  devices 
with  varying  channel  length  (defined  by  the  top  gate  width)  along  the 
same  bilayer  graphene  domain  wall.  Figure  4b  displays  the  measured 
conductance  of  the  two  sets  of  bilayer  graphene  domain  wall  channels 
with  lengths  of  200  nm,  400  nm,  600  nm  and  800  nm,  respectively  and 
another  batch  of  1,000-nm-long  bilayer  graphene  domain  wall  channels. 
All  the  conductance  values  are  of  the  order  of  the  quantum  conductance 
(2 e2/h,  where  e  is  the  unit  charge  and  h  is  Planck’s  constant),  and  one  of 
the  shortest  channels  exhibits  a  conductance  approaching  4 e2/h. 

The  observed  one-dimensional  conducting  channels  correspond  to 
the  chiral  edge  modes  at  AB-BA  domain  walls  recently  predicted  by 
theory11-13.  Figure  4c  shows  the  electronic  band  structure  around  the 
K  and  K'  points  for  bilayer  graphene  with  an  AB-BA  domain  wall.  In 
gapped  bilayer  graphene,  the  AB-  and  BA-stacked  domains  represent 
two  different  quantum  valley  Hall  insulators.  Consequently,  a  domain 
wall  separating  insulating  AB  and  BA  domains  will  confine  a  pair  of  one¬ 
dimensional  conducting  channels  along  the  boundary  for  each  valley. 
The  conducting  electrons  in  different  valleys  have  opposite  energy  dis¬ 
persions,  and  therefore  propagate  with  opposite  velocities.  Such  counter- 
propagating,  valley-projected  chiral  electrons  along  the  AB-BA  domain 
walls  are  illustrated  in  Fig.  4d  in  real  space.  These  conducting  chiral 
modes  are  topologically  protected  within  a  single  valley.  Inter-valley 
scatterings  can  potentially  mix  the  forward-  and  backward-moving  elec¬ 
trons  and  eliminate  the  chiral  modes,  analogous  to  magnetic  impurities 
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that  flip  spins  and  eliminate  chiral  modes  in  quantum  spin  Hall  insula¬ 
tors.  However,  at  atomic-defect-free  AB-BA  domain  walls  where  the 
crystal  lattice  is  smoothly  varying,  such  inter-valley  scattering  will  be 
absent  because  it  requires  a  large  momentum  change  of  4n /  (3\/3)  ac_  c, 
where  ac_c  is  the  carbon-carbon  bond  length  of  1.42  A. 

We  can  quantitatively  compare  the  observed  one- dimensional  con¬ 
ductance  along  the  AB-BA  domain  wall  with  theoretical  predictions. 
As  illustrated  in  Fig.  4d,  two  conducting  chiral  modes  are  present  at  the 
domain  wall  for  a  given  current  flow  direction.  In  addition,  each  chiral 
mode  is  doubly  degenerate  owing  to  electron  spin.  Therefore,  four  one¬ 
dimensional  conducting  channels  exist  at  the  domain  wall,  giving  an 
ideal  conductance  of  4 e2/h,  that  is,  two  conductance  quanta. 
Experimentally,  we  observed  conductance  of  the  order  of  the  conduc¬ 
tance  quanta  in  all  devices,  and  one  200-nm-long  domain  wall  channel 
exhibits  conductance  close  to  the  value  of  4e2/h  (Fig.  4b).  We  found 
that  overall  the  one-dimensional  channel  conductance  becomes  lower 
with  increasing  channel  length.  The  observed  length- dependent  con¬ 
ductance  through  bilayer  graphene  domain  walls  can  be  described  by 
the  Landauer-Biittiker  formula28  with  a  finite  electron  mean  free  path: 

G  =  (rc  +  Ro  ^1  +  ^  ,  where  L  is  the  channel  length,  R0  =  hi 4c2, 

the  contact  resistance  Rc  =  0,  and  the  mean  free  path  L0  =  420  nm 
(solid  line  in  Fig.  4b).  The  electron  mean  free  path  of  420  nm  can  arise 
from  a  finite  inter- valley  scattering  along  the  bilayer  graphene  domain 
walls,  and  previous  transport  studies  have  reported  a  comparable 
inter- valley  scattering  length  of  300-800  nm  in  bilayer  graphene  on 
Si02/Si  substrates29.  Nevertheless,  this  electron  mean  free  path  along 
the  bilayer  graphene  domain  wall  is  orders  of  magnitude  longer  than 
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Figure  4  |  Conductivity  of  topologically  protected  chiral  modes  at  AB-BA 
domain  walls,  a,  CNP  resistances  as  a  function  of  displacement  field  D  for 
bilayer  graphene  with  no  domain  wall  (NDW,  black  squares)  and  for  bilayer 
graphene  with  a  domain  wall  (DW,  red  circles)  obtained  from  Fig.  3.  The 
saturation  of  the  resistance  at  a  much  lower  value  at  high  D  originates  from  one¬ 
dimensional  conducting  channels  at  the  domain  wall,  b,  Length  dependence  of 
the  bilayer  graphene  domain  wall  conductance.  Green  and  blue  symbols 
represent  two  sets  of  devices,  each  fabricated  on  one  single  domain  wall  in  the 
same  bilayer  graphene  flake.  Red  symbols  represent  another  batch  of  1,000-nm- 
long  domain  wall  channels.  All  devices  exhibit  conductance  of  the  order  of  the 
quantum  conductance,  with  lower  conductance  in  longer  channels  and  a 
maximum  conductance  approaching  4 e2/h.  This  length- dependent  one¬ 


dimensional  conductance  can  be  described  by  four  conducting  channels  with  a 
mean  free  path  of  420  nm  using  the  Landauer-Biittiker  formula  (solid  line), 
c,  Electronic  band  structure  of  bilayer  graphene  with  an  AB-BA  domain  wall. 
Solid  green  and  orange  regions  represent  states  of  the  bulk  bilayer  with  a 
bandgap  of  A,  £F  indicates  the  Fermi  energy,  and  the  blue  (red)  lines  represent 
topologically  protected  K- valley  (K' -valley)  chiral  electron  modes  at  the 
domain  wall.  Each  chiral  electron  mode  is  doubly  spin  degenerate,  leading  to 
four  valley-polarized  one- dimensional  conducting  channels  at  the  AB-BA 
domain  wall  of  a  gapped  bilayer  graphene  sample,  d,  Illustration  of  valley- 
polarized  electrical  current  at  the  bilayer  graphene  AB-BA  domain  wall. 
Boundary  states  in  the  K  and  K'  valley  move  in  opposite  one-way  channels, 
leading  to  ballistic  charge  and  valley  current  along  the  domain  wall. 


that  in  the  bilayer  graphene  domain,  which  has  a  mobility  of 
—2,000  cm1 2  V- 1  s_1  and  a  mean  free  path  of  only  —8  nm  in  our  sam¬ 
ples.  In  general,  scattering  and  localization  effects  tend  to  be  more 
important  in  one  dimension.  The  much  longer  mean  free  path  at  the 
AB-BA  domain  walls  therefore  highlights  the  topological  protection  of 
chiral  electron  modes,  where  the  electron  velocity- valley  locking  for¬ 
bids  intra-valley  backscattering.  Because  inter- valley  scattering  along 
the  smooth  domain  walls  is  negligible,  the  chiral  valley  channels  at  AB¬ 
BA  domain  walls  are  very  robust  even  for  bilayer  graphene  devices  on 
Si02/Si  substrates.  We  envision  that  encapsulation  in  boron  nitride 
layers30  could  potentially  lead  to  chiral  electrons  at  AB-BA  domain 
walls  with  very  long  coherent  lengths. 

Our  observation  of  topologically  protected  one- dimensional  con¬ 
ducting  channels  at  bilayer  graphene  domain  walls  opens  up  exciting 
opportunities  to  explore  tunable  topological  phases  and  valley  physics 
in  graphene.  Owing  to  the  velocity- valley  locking,  the  one- dimensional 
chiral  modes  at  the  AB-BA  domain  walls  lead  not  only  to  ballistic 
charge  transport,  but  also  to  ballistic  valley  transport.  Such  chiral 
modes  can  enable  novel  ways  to  control  the  valley  degree  of  freedom 
for  valleytronic  devices  based  on  bilayer  graphene. 

Online  Content  Methods,  along  with  any  additional  Extended  Data  display  items 
and  Source  Data,  are  available  in  the  online  version  of  the  paper;  references  unique 
to  these  sections  appear  only  in  the  online  paper. 
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METHODS 

Mechanism  of  the  infrared  contrast  for  domain  walls.  Our  near -field  infrared 
nanoscopy  probes  the  local  optical  conductivity  at  infrared  frequencies  with  —40  nm 
spatial  resolution.  The  infrared  contrast  of  the  domain  wall  arises  from  its  different 
electronic  structure  and  infrared  conductivity  from  the  AB-stacked  bilayer  graphene. 
An  exact  description  of  the  local  band  structure  of  a  domain  wall  is  beyond  the  scope 
of  this  work.  However,  the  domain  wall  region  is  characterized  by  shifted  top  and 
bottom  layers,  which  can  be  approximated  by  two  randomly  stacked  monolayers  of 
graphene.  Such  randomly  stacked  graphene  has  electronic  structure  similar  to  that  of 
two  separated  monolayers  close  to  the  Dirac  point31,  and  its  infrared  conductivity  can 
be  approximated  by  2<7monoiayer.  The  AB-stacked  bilayer  graphene,  on  the  other  hand, 
has  an  infrared  conductivity  of  ObUayer-  The  optical  conductivities  of  both  monolayer 
and  AB-stacked  bilayer  graphene  have  been  well  established  in  the  literature32.  At  the 
infrared  wavelength  of  6.1  pm,  the  value  of  2cmonolayer  is  about  20%  higher  than  that 
of  g biiayen  which  can  qualitatively  explain  the  bright  infrared  contrast  of  the  domain 
walls  observed  in  our  studies. 

Electrical  transport  data  of  four  extra  bilayer  graphene  devices.  We  investi¬ 
gated  five  bilayer  graphene  device  pairs  with  a  fixed  channel  length  of  1  pm.  In 
Extended  Data  Fig.  1,  we  show  the  complete  set  of  data  taken  from  the  four  pairs  of 
bilayer  graphene  devices  other  than  that  included  in  Fig.  3c  and  d.  Each  pair  of 
bilayer  graphene  devices  includes  one  with  AB-BA  stacking  domain  walls  and  one 
reference  device  with  no  domain  walls  on  the  same  graphene  flake.  Extended  Data 
Fig.  la-d  displays  resistances  of  corresponding  bilayer  graphene  devices  with 
domain  walls,  including  three  single-domain-wall  devices  (Extended  Data  Fig. 
la-c)  and  one  double- domain- wall  device  (Extended  Data  Fig.  Id).  All  single¬ 
domain-wall  devices  show  resistances  saturating  at  14  kQ  with  increasing  D,  and 
the  double- domain-wall  device  has  a  resistance  saturating  around  8  kQ.  For  com¬ 
parison,  Extended  Data  Fig.  le-h  displays  the  gate-dependent  resistance  for  four 
bilayer  graphene  devices  without  boundaries.  These  devices  show  continuously 
increasing  CNP  resistance  to  well  above  14  kQ  (8  kQ  in  Extended  Data  Fig.  lh),  as 


indicated  by  the  dashed  line  at  high  Vbg.  All  the  data  show  that  one- dimensional 
conducting  channels  are  present  at  AB-BA  domain  walls,  and  that  they  dominate 
electrical  transport  in  gapped  bilayer  graphene. 

Length  dependence  of  electrical  transport  along  AB-BA  domain  walls.  We 

investigated  two  sets  of  dual-gated  bilayer  graphene  domain  wall  devices  of  vary¬ 
ing  channel  lengths  (defined  by  the  top  gate  width).  Each  set  of  devices  are 
fabricated  on  one  single  domain  wall  in  the  same  bilayer  graphene  flake,  and 
includes  four  channels  with  lengths  of  200  nm,  400  nm,  600  nm  and  800  nm, 
respectively.  The  complete  transport  data  of  these  bilayer  graphene  domain  wall 
devices  are  shown  in  Extended  Data  Fig.  2. 

Temperature  dependence  of  electrical  transport  of  AB-BA  domain  walls.  We 

investigated  the  temperature  dependence  of  electrical  transport  both  through 
AB-stacked  bilayer  graphene  and  through  AB-BA  domain  walls.  Extended 
Data  Fig.  3  displays  representative  data  for  CNP  resistance  at  Vbg  =  —140V 
for  a  bilayer  graphene  device  without  a  domain  wall  (black  line)  and  a  device 
with  a  domain  wall  (red  line).  For  the  bilayer  graphene  bulk  device  without  a 
domain  wall,  the  CNP  resistance  increases  by  over  300  kD  (2.5  times)  as  the 
temperature  is  varied  from  50  K  to  1.8  K.  This  large  increase  of  resistance  at 
lower  temperatures  is  characteristic  of  insulating  behaviour.  In  comparison, 
the  CNP  resistance  of  device  with  a  domain  wall  (red  line)  has  much  weaker 
temperature  dependence:  the  resistance  increase  is  only  5kQ  (—50%),  and 
part  of  the  increase  is  due  to  the  contribution  from  bulk  bilayer  graphene 
channel  in  parallel  with  the  one-dimensional  domain  wall  channel.  The  weak 
temperature  dependence  of  conductance  along  the  domain  wall  is  character¬ 
istic  of  a  metallic  system. 

31.  Lopes  dos  Santos,  J.  M.  B.,  Peres,  N.  M.  R.  &  Castro  Neto,  A.  H.  Graphene  bilayer 
with  a  twist:  electronic  structure.  Phys.  Rev.  Lett.  99,  256802  (2007). 

32.  Mak,  K.  F.,  Sfeir,  M.  Y.,  Misewich,  J.  A.  &  Heinz,  T.  F.  The  evolution  of  electronic 
structure  in  few-layer  graphene  revealed  by  optical  spectroscopy.  Proc.  Natl 
Acad.  Sci.  USA  107, 14999-15004  (2010). 


©2015  Macmillan  Publishers  Limited.  All  rights  reserved 


RESEARCH 


LETTER 


Extended  Data  Figure  1  |  Electrical  transport  of  dual-gated  bilayer 
graphene  devices  with  (left)  and  without  (right)  AB-BA  domain  walls. 

a-d,  Gate-dependent  resistance  R  for  three  bilayer  graphene  devices  with  one 
domain  wall  (a-c)  and  one  with  two  domain  walls  (d)  at  4.2  K.  The  backgate 
voltage  Vbg  is  varied  from  0  V  to  140  V  with  a  step  of  10  V,  and  the  top  gate 
voltage  Vtg  is  swept  continuously.  The  resistance  peak  in  each  trace  corresponds 
to  a  CNP.  The  resistance  at  the  CNP  first  increases  and  then  saturates  at 


~14kQ  in  a-c  (~8  kQ  in  d).  The  dashed  lines  indicate  the  saturated 
conductance  values,  e-h,  Gate-dependent  resistance  of  four  reference  bilayer 
graphene  devices  without  domain  walls.  The  resistance  at  the  CNP  increases 
continuously  at  high  Vbg  owing  to  the  opening  of  a  bandgap,  and  reaches  values 
much  higher  than  that  of  corresponding  devices  with  bilayer  graphene  domain 
walls  (indicated  by  the  dashed  line). 
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Extended  Data  Figure  2  |  Length  dependence  of  electrical  transport  of  AB¬ 
BA  domain  walls.  Data  in  left  panels  are  from  one  graphene  flake,  data  in  right 
panels  are  from  another,  a-d,  Gate- dependent  resistance  for  four  bilayer 
graphene  devices  fabricated  on  a  long  domain  wall  with  channel  lengths  of 
200  nm,  400  nm,  600  nm  and  800  nm  (defined  by  the  top  gate  width), 
respectively,  e-h,  Gate-dependent  resistance  for  four  bilayer  graphene  devices 
fabricated  on  another  long  domain  wall  with  channel  lengths  of  200  nm, 


400  nm,  600  nm  and  800  nm  respectively.  The  backgate  voltage  Vbg  is  varied 
from  0  to  —140  V  with  a  step  of  10  V,  and  the  top  gate  voltage  Vtg  is  swept 
continuously.  The  resistance  peak  in  each  trace  corresponds  to  a  CNP.  The 
resistance  at  the  CNP  first  increases  as  the  bandgap  is  opened  in  the  bulk  and 
then  saturates  at  high  Vbg,  where  the  electrical  transport  is  dominated  by  one¬ 
dimensional  channels  along  the  domain  wall. 
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